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Abstract 

An analysis of the acoustical and flowfield environment for the scaled 1 -pound-force 
(Ibf) thrust tabletop motor was performed. The jet characterization is based on computa- 
tional fluid dynamics (CIT)) in conjunction with Kirchhoff surface integral formulation and 
compared with correlations developed for measured rocket noise and a pressure fluctuation 
scaling (PFS) method. Comparisons are made for the overall sound pressure levels 
(OASPL’s) and spectral dependence of sound pressure level (SPL). 

INTRODUCTION 

At NASA Kenned) Space Center, a Launch Systems Testbed (LST) is currently under 
development to establish a capability to simulate small-scale launch vehicle environment for 
use in testing and evaluation of launch pad designs for future space vehicles. These predic- 
tion methods include computational fluid dynamics (CFD), analytical correlations, and pres- 
sure fluctuation scaling (PFS) methods. This report, primarily derived from [1], summarizes 
the analytical studies of a small-scale tabletop rocket that is being acquired from NASA 
Stennis Space Center for the purpose of establishing the initial test facility and instrumenta- 
tion. 
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DESCRIPTION OF THE TABLETOP ROCKET 

Fig. 1 shows a schematic of the nozzle geometry tor the tabletop motor [2] . This motor 
has cylindrical Plexiglas fuel burning in gaseous oxygen. The combustion gas is approxi- 
mated as carbon dioxide (CO 2 ). The nozzle exit conditions, as derived form isentropic ex- 
pansion [3] are shown in Table 1. 


ANALYSIS 


CFD/Kirchhoff Analysis 

The CFD/Kirchhoff analysis is based on the application of OVERFLOW CFD Navier- 
Stokes code [4, 5] for identifying the noise sources in the nonlinear source field and 
Kirchhoff surface integral [6] for the propagation of sound radiation to the near field and the 
far field. For the acoustic radiation, the Kirchhoff code YORICK [7] was considered. 

CFD Solution. For the CFD computations, an axisymmetric grid of 200 x 100 is considered. 
The CFD solution converged after about 15,000 time-step iterations (the code is run in a 
time-accurate manner) before a periodic state is established. 

Acoustic Solution. After a periodic state is established, the data from CFD is communicated 
to the Kirchhoff code. The radius of the Kirchhoff surface is taken 6 radii from the jet axis. 

Correlation/Analytical Method 

Leneman [8] proposed an empirical correlation for the OASPL as 

OASPL(dB) = 115 + 10 log(F/)-201og(x) (la) 

where F denotes the total thrust (Ibf), / the specific impulse (sec), and x the axial distance 
(ft) in the far field from he jet exit plane. Margasahayam and Caimi [9] proposed as: 

OASPL (dB) = 1 1 5 + 1 0 log(F /) - 20 log(jc) - 20x ( 1 b) 

with the fourth term providing a near-field correction to the Leneman correlation. 

The NASA SP-8072 method [10] represents a more accurate correlation compared to 
the Leneman correlation and is based on normalized relative power spectrum from test data. 
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Pressure Fluctuation Scaling Method 

In the method of pressure fluctuation scaling due to Norton et al. [11], which has been 
successfully used for industrial gas piping networks, etc., the nondimensional fluctuating 
pressure spectral density d>p is given by 

4G„(ty) 

O ra) - f \ , Q. = (Miu (2) 

p^U^a 

where Q. is the nondimensional frequency (Strouhal number), Gp is the fluctuating pressure 
spectral density, p is the gas density, U is the flow velocity, a is the nozzle exit radius, and 
CO = iTtf. The turbulent wall pressure spectrum (shell or flat plate) provides a lower level, 
and the severe internal flow disturbance provides the upper level for scaling purposes. 

The sound spectrum is then calculated at the exit downstream of the shock. Power 
spectral estimates of fluctuating pressure for the tabletop rocket are obtained from 




^(Q.)p'^U^a^co, 


(D: - UQ.:la 


(3) 


where AttJi is the nominal bandwidth of band i with center frequency cOi. Coirection factors 
are then applied for the o])en jet exhausting out of the nozzle. At this stage, a 1/r (r is the dis- 
tance from the sound source) decay with distance from the nozzle is assumed for circumfer- 
ential radiation away from the exit plane (e.g., position A, Fig. 1). The 1/r decay is also con- 
sidered for post shock centerline decay (e.g., position X) because the mean flow will en- 
hance/convect the monojDole- and dipole-type radiation from the duct exit. These assump- 
tions need to be refined. 


RESULTS AND COMPARISON 


Flowfleld Results 

The axial variatior of Mach number at the centerline (Fig. 2) indicates that the core 
length of the jet extends to about 25 jet radii from the exit plane. The average Mach number 
in the core region downstream of the shock is about 0.6, characterizing the jet as subsonic. 

Acoustical Results 

Table 2 summarize s the comparison of OASPL at stations A, E, and X (Fig. 1) as esti- 
mated by the different methods. 
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The spectral composition of 1/3 octave SPL at point A in the exit plane is presented in 
Fig. 3a. The SPL from the upper-level PFS exceeds the CFD result. However, both the CFD 
and the upper-level PFS method suggest a peak frequency of about 20 kilohertz (kHz), which 
corresponds to a Strouhal number St = 0.20. The lower-level PFS method produces SPL that 
is close to that of the CFD but produces a peak frequency of about 100 kHz, which appears 
somewhat overestimated. 

A comparison of the 1/3 octave SPL content at location E on the axis downstream of 
the exit plane is portrayed in Fig. 3b. In general, the SPL values at E are considerably higher 
than those at position A, due to the near-field location of E and the effect of directivity. A 
peak frequency of aboul 35 kHz (corresponding to St = 0.35) is indicated by CFT), the 
NAS A-SP method, and the PFS upper-level method. The CFD results for the SPL are higher 
than those from the NAS A-SP method. In the case of the PFS method, the upper-level val- 
ues for the SPL are closer to the CFD, whereas the lower-level results match well with the 
NASA-SP method. 

Fig. 3c compares the spectral distribution of 1/3 octave SPL at location X, which is 
downstream of location E but 1 jet radius away from the axis. In general, the sound pressure 
levels at X are smaller than those at E. A reduced peak frequency of about 25 kHz (corre- 
sponding to St = 0.25), relative to location E, is indicated by CFD, the NASA-SP method, 
and the PFS upper-level method. The SPL variation between the various methods at position 
X is somewhat smaller compared to that at positions A and E. 

A comparison of the 1/3 octave SPL at locations A, E, and X as predicted by CFD is 
shown in Fig. 4. Downstream of the exit plane, the peak frequency level is increased in the 
near field, as the distance from the exit plane is decreased. While there appears to be a simi- 
larity in the shape of the ,SPL spectrum between E and X, there is a departure in similarity at 
position A for frequenci(JS in excess of about 40 kHz, where a rapid decline in SPL is ob- 
served. This lack of similarity in the SPL spectrum at position A is perhaps attributable to 
the existence of a norma shock near the nozzle exit or to the far-field directivity effects. 

CONCLUSIONS 

The CFD/Kirchhoff method provided flowfield structure and distributions of the 
OASPL and SPL spectmim for the tabletop rocket. Because of the very small size of the noz- 
zle exit diameter, the jet is distinguished by high peak frequency levels (20 to 35 kHz) for the 
SPL spectrum. A comparison of the results from CFD, the NASA-SP method, and the PFS 
method at three spatial locations (A, E, and X) showed that a similarity in shape is generally 
noted for the SPL composition. There is some disagreement between the different methods 
for the peak frequencies and the SPL values. A more realistic assessment of these methods 
can only be carried out w ith the aid of an experimental test program currently under devel- 
opment. 
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Table 1 . Summary of Nozzle Parameters 


Parameter 

Value 

Stagnation pressure, psia 

54.7 

Stagnatitm temperature, °R 

1960 

Nozzle mass flow rate, Ibm/s 

0.0247 

Nozzle throat diameter, inch 

0.2 

Exit pressure, psia 

2.38 

Exit temperature, °R 

987 

Exit velocity, ft/s 

3,170 

Nozzle exit Mach number 

2.65 

T 

Nozzle exit density, Ibm/ft 

0.0099 

Jet exit Reynolds number 

2.9x10"^ 


Table 2. Comparison of QASPL at Locations A, E, and X 


Location 

CFD 

Leneman 


PFS 


Modified 

Upper Level 

Lower Level 

A 

118 


— 

— 

153 

119 

E 


175 

155 

125 

176 


X 

140 

144 

124 

117* 

135 



* Not applicable, ** Value corresponding to Z = 0 



Fig. 1. Schematic of the Tabletop Rocket Nozzle 
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Fig. 2. Jel Centerline Variation of Mach Number From CFD 



frequency (kHz) 


Fig. 3a. Comparison of 1/3 Octave Spectra at Location A 
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Fig. 3c. Comparison of 1/3 Octave SPL Spectra at Lx)cation X 



Fig. 4. Comparison of 1/3 Octave SPL Spectra From CFD at A, E, and X 






